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Summary. Complementary oligonucleotide primers which flank a 1146-nucleo¬ 
tide gene fragment (SIB: nt 1185 to 2333) encompassing a polymorphic region 
(nt 1368 to 1776) of the SI subunit of bovine coronavirus spike glycoprotein 
were used for enzymatic amplification by PCR. We chose four clinical isolates, 
recovered from cases of epidemic diarrhea in neonatal calves in Quebec dairy 
herds between 1987-1990, to specifically amplify and analyze their sequences 
in the selected genomic area. Nucleotide sequence analysis of the four clinical 
isolates indicated that their SIB gene fragments were highly conserved. We 
also compared the SIB gene sequences of the Quebec BCV isolates to the 
published corresponding sequences from BCV-L9 [37], BCV-MEB [1], and 
BCV-F15 [3] reference strains. A high degree of similarity was demonstrated 
for all viruses, no deletions or insertions were observed, and the only variations 
that were identified consisted of nucleotide substitutions. The differing nucleotides 
and amino acids (aa) were not distributed randomly over the entire sequence 
but rather were clustered in the polymorphic region. Of these, four sporadic 
aa changes were located in antigenic domain II (aa residues 517 to 720) of SI. 
This correlates with varied antigenicity observed among the BCV Quebec 
isolates when reacting with MAbs directed against the S glycoprotein of the 
Mebus strain. The other mutations seem to be fixed in all Quebec isolates. 

Introduction 

Bovine coronavirus (BCV) is an important cause of acute enteritis in newborn 
calves [18]. Some recent reports suggest that BCV is also responsible for upper 
respiratory tract illness in growing calves [22,24], winter dysentery (hemorrhagic 
enteritis) or chronic shedding in adult cattle [2, 25]. The viral particle is mostly 
spherical, enveloped with a diameter of around 100 nm, and displays two fringes 
of surface projections [28, 31]. The viral genome consists of a large single- 
stranded RNA with positive polarity, approximately 30 Kb in length [31], and 
encodes four major structural proteins which are the phosphonucleoprotein 
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(N; 52kDa), the matrix glycoprotein (M; 24-26 kDa), the spike protein (S; 
180-200 kDa) and the hemagglutinin/esterase (HE; 62-65 kDa) [11, 13, 16]. 

The S glycoprotein often is post-translationally cleaved by host-cell pro¬ 
teases into two 90-100kDa fragments, SI and S2, respectively corresponding 
to the N- and C-terminal subunits [3]. This proteolytic cleavage may have an 
enhancing effect on viral infectivity and fusion activity [5, 29, 30]. Both S 
subunits contain antigenic domains responsible for inducing neutralizing 
antibodies [9, 10, 15, 33], but SI appears to elicit the production of the 
monoclonal antibodies (MAbs) displaying the highest neutralizing activity [34]. 
Competitive ELISA tests using anti-S MAbs suggest the existence of at least 
four independent antigenic domains (A, B, C and D) for the SI subunit, antigenic 
sites S1A and SIB being associated to the neutralizing activity [34], Alterations 
in the S protein resulted in changes in viral antigenicity, as demonstrated for 
neutralization-escape mutants selected under the pressure of MAbs specific for 
the S protein of BCV [36]. A highly polymorphic region of the SI subunit of 
BCV (aa residues 456 to 592) and the mouse hepatitis virus type 4 (MHV-4) 
(aa residues 426 to 477) has been identified [20, 21, 36]. This region has been 
demonstrated to undergo deletions or point mutations selectable by neutralizing 
MAbs [21]. The mutant viruses displayed decreased virulence, suggesting that 
the SI polymorphic region may be directly involved in the viral pathogenicity 
[ 12 ]. 

Recent comparison of Quebec BCV isolates using polyclonal antisera and 
MAbs directed against the S glycoprotein of the Mebus strain confirmed their 
close antigenic relationship, but also revealed that they are assigned to at least 
three distinct antigenic subgroups [19]. In the present study, we analysed the 
nucleotide (nt) sequence of a 1146-nt gene fragment (SIB) encompassing the 
highly polymorphic region and the proteolytic cleavage site of 4 serologically 
distinct Quebec BCV isolates and compared them to the homologous sequences 
of the Mebus [1] and FI5 [3] reference strains. 

Materials and methods 

Cells and viruses 

The cell culture-adapted Mebus strain [18] of BCV was obtained from the American Type 
Culture Collection (ATCC VR-874), Rockville, MD. This prototype BCV strain was 
originally isolated in bovine fetal kidney cells from diarrhea fluid of calf [18]. Viral isolates 
BCQ9, BCQ 571, BCQ 20 and BCQ 2070 were recovered during winter 1989 from clinical 
cases of epidemic diarrhea in newborn calves affecting Quebec dairy herds located in four 
different geographic areas [19]. No commercial BCV vaccine had been applied in these 
herds during the year preceeding the emergence of clinical cases. The various Quebec BCV 
isolates could be differentiated serologically by indirect immunofluorescence, hemagglutina¬ 
tion-inhibition, seroneutralization and competitive ELISA tests using Mabs directed against 
the S protein of the prototype Mebus strain [19]. For the present studies, the BCV strains 
were passaged not more than five times in human rectal tumor (HRT-18) cells in the 
presence of lOU/ml of bovine pancreatic trypsin [8]. The extracellular virons were purified 
from the supernatants of infected cell cultures by differential and isopycnic ultracentrifugation 
on sucrose gradients as previously described [8], 
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Isolation of viral RNA 

Viral RNA was extracted following the guanidinium isothiocyanate-acid phenol method 
[4]. Briefly, 100 pi of purified virus was resuspended and vortexed vigorously in 500 pi of 
solution D (4M guanidinium isothiocyanate, 0.75 M sodium citrate, 0.5% Sarkosyl and 
0.1 M 2-mercaptoethanoi). After acid phenol treatment, the aqueous phase was removed 
and RNA was precipitated with isopropanol. The RNA pellets were washed with 75% 
ethanol, dissolved in 20 pi RNase-free water, and stored at — 70 °C. 


Preparation of oligonucleotide primers 

The oligonucleotide primers used in the RT-PCR assay were chose on the basis of 
the previously published sequence for the Mebus strain of BCV [1], EMBL/Gen Bank 
accession No. M 31053. The primers were: oligo Sis (5T GCT AT ACCG A AT GGT AGG 3 ? ) 
complementary to a sequence at positions 1185-1203 on the S gene and oligo Si as 
(5 GTAAACCGATAACCAGTGG3') corresponding to antisense of the sequence at position 
2315-2333 on the S gene. This amplification primer pair flanks a 1146-nucleotide gene 
fragment (SIB) encoding most neutralizing epitopes and the proteolytic cleavage site [34,36]. 
In order to clone the PCR products using the CloneAmp system (GIBCO BRL, Gaithersburg, 
MD), the following sequences were incorporated into the 5' end of each amplification 
primer: Sn 5-CUA CUA CUA CUA-3', and Asn 5'-CAU CAU CAU CAU-3'. Oligonucleo¬ 
tide syntheses and incorporation of dTMP into 5' termini were carried out on an automated 
DNA synthesizer (Pharmacia Biotech Inc., Baie d’Urfe, QC, Canada). 


First-strand cDNA synthesis 

Prior to PCR amplification, the target RNA was converted to a cDNA by use of reverse 
transcriptase (RT) and SI as primer. The RT reaction contained 2 pi of 10X PCR buffer 
(500mM KC1, lOOmM Tris-HCl (pH 9), 15mM MgCl 2 , and 1% Triton X-100), lpl 
containing 20 units of RNA guard (Pharmacia Biotech Inc.), 2 pi of 10 mM each dNTP, 
2 pi of 25 mM MgCl 2 , 1 pi of 50pmol antisense primer, and 2-5 pi of the extracted viral 
RNA. A 20-pl total reaction volume was obtained by adding sterile DEPC-treated distilled 
water. After heating to 94 °C for 1 min and 65 °C for 5 min, 20 units MuLV RT were added 
(Pharmacia Biotech Inc.). The mixture was incubated 60 min at 42 °C, and the RT reaction 
was stopped by heating to 95 °C for 5 min and stored at —20°C. 


PCR amplification of first-strand cDNA 

A volume of 80 pi of the PCR reaction mix was added to 20 pi of the cDNA mix such that 
the final 100pi volume contained lOmM Tris-HCl (pH 8.3), 50mM KC1, 2.5mM MgCl 2 , 
0.01% gelatin, 50pmol of Sis and Slas primers and 2.5 units of Taq DNA polymerase 
(BIO 101, La Jolla, CA, U.S.A.). After overlaying with 100 pi of mineral oil, amplification 
was carried out in a DNA thermal cycler (Perkin Elmer) with the following cycling program: 
94 °C denaturation for 1 min, 60 °C annealing for 1 min, and 72 °C elongation for 2 min. 
A total of 35 cycles was performed, the final polymerisation step was at 74 °C for 10 min. 
Following amplification, 5 pi of the final product was analyzed on a 0.8% agarose gel in 
TBE buffer (0.089 M Tris-HCl, (pH 8.3), 0.089 M boric acid, and 2 mM EDTA) and stained 
with ethidium bromide. 
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Cloning of PCR-generated fragments 

The PCR products were ligated into the pAMP 1 vector using the CloneAmp system 
(GIBCO BRL). The ligation mixture was used to transform competent cells E. coli DH5a 
(GIBCO BRL) following the manufacturer’s direction. 

Sequencing of cloned PCR fragment 

DNA sequencing was performed by the dideoxynucleotide chain-termination method [26], 
following the procedure outlined in the instruction of T7 sequencing kit (Pharmacia Biotech 
Inc.). Sequences of the additional primers used for sequencing the SIB fragment were: 5'- 
CGTGTAAATTGGATG-3', 5'-GTTGGCAAGTACAAG-3', 5'-AATCCTTGTACTTGC- 
3', 5'-GTTAAGTAGTCTCTA-3', 5'-AATTTACACGGACAG-3' and 5-TAGAGACTA- 
CTTAAC-3'. To assess the error rate of the reverse transcriptase and TaqDNA polymerase, 
clones from different cDNA synthesis and PCR events were sequenced. Errors generated 
by T7 polymerase were avoided by sequencing clones in both DNA strands and directions. 
Sequences were analysed by the computer programs, McVector 3.5 (Rainbow Tech. Inc., 
CA, U.S.A.) and GeneWorks 2.2 (IntelliGenetics Inc., Mountain View, CA, U.S.A.). 

Nucleotide sequence accession numbers 

The nucleotide sequence accession numbers (EMBL/GeneBank/DDBJ) are as follows: 
U06090 for BCQ-2070; U06093 for BCQ-571; U06092 for BCQ-20; and U06091 for BCQ-9. 

Results 

Comparison of SIB sequences of various clinical BCV isolates 

The complete SIB nucleotide (nt) and deduced amino acid (aa) sequences of 
the four clinical BCV isolates were aligned and compared (Figs. 1 and 2). Each 
SIB gene fragment is 1146bp long encoding the antigenic domain II identified 
within aa residues 517 to 720 [36], the polymorphic region located between 
aa residues 456 and 592 [36], and the proteolytic cleavage site beginning with 
aa 763 [1, 36]. The degree of homology and sequence conservation in the 
nucleotide and deduced aa sequences of the SIB gene fragment were analysed 
using a microcomputer program GeneWorks. The isolates BCQ-2070, BCQ-9 
and BCQ-571 possessed at least 98% nt sequence identity in paired comparisons. 
However, nt sequence of the SIB gene fragment of BCQ-20 was only 96% 
conserved and, thus, was the most distantly related. These sequences demonstrated 
a high degree of similarity, and the only variations observed consisted of point 
mutations that apparently not occurred at random. Fifty percent of the 
nucleotide substitutions were located in the third positions of the respective 
codons and therefore represent translationally silent mutations. This resulted 
in only thirteen aa changes among the four isolates. The most similar deduced 
aa sequences were those of BCV-2070 and BCV-571 or BCV-2070 and BCV-9 
(5 aa differences) in paired comparisons, and the most distant were those of 
BCV-571 and BCV-9 (10 aa differences). Four potential N-linked glycosylation 
sites were predicted in cases of the SIB fragment of the four Quebec BCV 
isolates, these glycosylation sites being located downstream the polymorphic 
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BCV-MEB TTTGCTATAC CCAATGGTAG GAAGGTTGAC CTACAATTGG GCAATTTGGG CTATTTGCAG TCTTTTAACT ATAGAATTGA TACTACTCCT ACAAGTTGTC 
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AGTTGTATTA TAATTTACCT GCTGCTAATG TTTCTGTTAG CAGGTTTAAT CCTTCTACTT GGAATAGGAG ATTTGGTTTT ACAGAACAAT TTGTTTTTAA 
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BCV-MEB OCCTCAACCT OTAOOTGTTT TTACTCATCA TGATGTTGTT TATGCACAAC ATTGTTTTAA AGCTCCCTCA AATTTCTGTC CGTGTAAATT GGATGGGTCT 300 

BCV-P15 .O.,.,..A.. . .. 300 

BCQ-2070 . ..C.G. *. *_._ *.A.*.. ... 300 

BCQ-571 .. C........G.... ......A.. .......... .......... 300 

BCQ-20 ... ,C. 4...aC*....... .......A.. .......... .......... 300 

BCQ-9 . ,C.G.... .A.... .A 300 


BCV-MEB 

BCV-F15 

TTGTGTGTAO GTAATGQTCC 
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BCV-MEB CTGTTCGGGT CTTGCTATTA AAAGTGATTA TTGTGGAGGT AATCCTTGTA CTTGCCAACC ACAAGCATTT TTGGGCTGGT CTGTTGACTC TTGTTTACAA 


BCV-F15 .... ...... I.*.. . 

BCQ-2070 .....C. .. ........C. .......... ... .....T...* 

BCQ-571 .......... ....... ...T.. . ..C 

BCQ-20 . ..C.G. 

BCQ-9 .T.. . <.C...T. 
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BCV-MEB GGGGATAGGT GTAATATTTT TGCTAATTTT ATTTTGCATG ATGTTAATAG TGGTACTACT TCTTCTACTG ATTTACAAAA ATCAAACACA GACATAATTC 700 

BCV-F15 . 700 

BCQ-2070 ...... . . ... . . . ... 4 *.,*,. 700 

BCQ-571 .*.*. 700 

BCQ-20 .C . 700 

BCQ-9 . 700 

BCV-MEB TTGGTGTTTG TGTTAATTAT GATCTTTATG GTATTACAGG CCAAGGTATT TTTGTTGAGG TTAATGCGAC TTATTATAAT AGTTGGCAGA ACCTTTTATA 800 

BCV-F15 -*..T.C.* . 800 

BCQ-2070 . .. . ..T* , . . ..... .. ,.,...,., 4 800 

BCQ-571 *.*..... 800 

BCQ-20 .T. 800 

BCQ-9 .T.. 800 

BCV-MEB TGATTCTAAT GGTAATCTCT ATGGTTTTAG AGACTACTTA ACAAACAGAA CTTTTATGAT TCGTAGTTGC TATAGCGGTC GTGTTTCAGC GGCCTTTCAT 900 

BCV-FI 5 ...... . ... ... . . . . . . 900 

BCQ-2070 .. .T.. ..... 900 

BCQ-S71 .!.C. 900 

BCQ-20 .T.\y.T. 900 

BCQ-9 ..* .T.T. 900 

BCV-MEB GCTAATTCTT CCGAACCAGC ATTGCTATTT CGGAATATTA AATGCAATTA CGTTTTTAAT AATACTCTTT CACGACAGCT GCAACCTATT AACTATTTTG 1000 

BCV-F15 ......4.., .. ,C...* ...... 1000 

BCQ-2070 .....C ........... 1000 

BCQ-571 C. 1000 

BCQ-20 .C.GC . . . *.. 1000 

BCQ-9 .C.GC.T.. 1000 

BCV-MEB ATAGTTATCT TGGTTGTGTT GTCAATGCTG ATAATAGTAC TTCTAGTGTT GTTCAAACAT GTGATCTCAC AGTAGGTAGT GGTTACTGTG TGGATTACTC 1100 

BCV-F15 4...C..4 .....G.C.. . 44, __ ... .. .. U00 

BCQ-2070 -C... 1100 

BCQ-571 ...T. 1100 

BCQ-20 ,...C . 1100 

BCQ-9 . . . .C..... 1100 

BCV-MEB TACA AAAAGA CGAAGTCGTA GA GCGATTAC CACTGGTTAT CGGTTTA 1147 

BCV-F15 ...CGTA. 1147 

BCQ-2070 T. 1147 

BCQ-571 . 1147 

BCQ-20 .GT. 1147 

BCQ-9 T. 1147 


Fig. 1. Nucleotide sequence comparisons of the SIB gene fragment (nt 1185 to 2333) of 
four Quebec BCV isolates, BCQ-2070, BCQ-571, BCQ-20 and BCQ-9, and two reference 
strains, BCV-MEB [1] and BCV-F15 [3]. Bold letters show the polymorphic region 
corresponding to aa residues 456 to 592. The underlined sequences refer to the proteolytic 
cleavage site. The sequence of reference strain BCV-MEB is written in full, while only 
changes in nucleotides are indicated for the sequences of the other strains 
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BCV-MEB GDRCNIFANF ILHDVNSGTT CSTDLQKSNT DIILGVCVNY DLYGITGQGI FVEVNATYYN SWQNLLYDSN GNLYGFRDYL TNRTFMIRSC YSGRVSAAFH 700 


BCV-F15 .A. 700 

BCQ-2070 700 

BCQ-571 .L. 700 

BCQ-20 . 700 

BCQ-9 . 700 

BCV-MEB ANSSE PALLF RNIKCNYV FM MT LSRQLQPI NYFDSYLGCV VNADNSTSSV VQTCDLTVGS GYCVDYS TKR RSRRA ITTGY RF 782 

BCV-F15 . , . .. “77 .A.A.V. 782 

BCQ-2070 .S. 782 

BCQ-571 .V. 782 

BCQ-20 .S.S. 782 

BCQ-9 .S.I.S. 782 


Fig, 2. Alignment of the deduced amino acid sequences of the SIB gene fragment of four 
Quebec BCV isolates, BCQ-2070, BCQ-571, BCQ-20, and BCQ-9; and two reference strains, 
BCV-MEB [1], and BCV-F15 [3], The polymorphic region (aa residues 456 to 592) is 
written in bold face, the proteolytic cleavage site (aa residues 663 to 767) is underlined and 
the potential N-linked glycosylation sites are double underlined. The sequence of reference 
strain BCV-MEB is written in full, while only changes in amino acid residues are indicated 

for the sequences of the other strains 


region and antigenic domain II. Twenty-six cysteine residues were conserved 
among all of the Quebec isolates. 

Comparison of SIB sequences of Quebec BCV isolates with 

other known BCV strains 

As illustrated in Fig. 1, the complete SIB nt sequences of the four Quebec 
clinical BCV isolates were aligned and compared with the homologous 
sequences of the Mebus (BCV-MEB) and F15 (BCV-F15) reference strains. The 
level of nt sequence identity between the SIB fragments of the four Quebec 
BCV isolates and that from BCV-MEB and BCV-F15 strains was 97% and 
96%, respectively. In contrast, the polymorphic region (nt 168-576) seemed to 
be more conserved among the Quebec BCV isolates and BCV-F15 (97%). 

Frameshift, deletion, or insertion, and non sense mutations were not ob¬ 
served; the only variations among these sequences consisted of 40 nt substitu¬ 
tions, which represented 4% of the sequence. The differing nucleotides were not 
distributed randomly over the entire sequence but rather were clustered in 
the polymorphic region. This was also reflected in the deduced aa sequences 
(Fig. 2). The 60 aa residues immediately upstream the polymorphic region were 
almost completely conserved among all the isolates. However, 12 aa substitutions 
were found between residues 456 and 592 of the polymorphic region. The fact, 
that these aa changes occurred in almost identical location suggests these 
differences may be significant. The 29 aa residues (aa 593 to 621) immediately 
downstream the hypervariable region have been shown to be part of the 
antigenic domain II which is recognized by neutralizing Mabs BB7-14 [36]. 
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This region was totally conserved among all BCV isolates tested. The twenty-six 
cysteine residues identified above were perfectly conserved among all of the 
Quebec isolates. The five remaining differences in the 160 residues upstream 
the proteolytic cleavage site represented variability which seemed to be unique 
to some isolates. The residue changes consisted of Phe for Leu and Asp for 
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Fig. 3. Hydrophilicity profile of the SIB gene fragment of four Quebec BCV isolates 
(BCQ-2070, BCQ-571, BCQ-20, and BCQ-9); and two BCV reference strains, BCV-MEB 
[1], and BCV-F15 [3]. Hydrophilicity plots (according to the Kyte-Doolittle methods, 
using a seven-amino acid window) starts at amino acid 401 of the mature SI protein. 
Arrows in hydrophilicity plots of isolates BCQ-20 and BCQ-9 show regions of increased 

hydrophilicity 
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8CV-VACC 
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BCV-MEB 


BCV-LY 
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Fig. 4. Phylogenetic tree of the SIB gene fragment of four Quebec BCV isolates (BCQ-2070, 
BCQ-571, BCQ-20, and BCQ-9); and six BCV strains; BCV-MEB [1]; BCV-F15 [3]; 
BCV-L9, BCV-LY, and BCV-VACC [37]; and BCV-QUEB [20], The length of the 
horizontal lines connecting one sequence to another is proportional to the estimated genetic 

distance between the sequences 


Val for BCQ-571, He for Ser for BCQ-20 and BCQ-9, and Thr for lie for 
BCQ-9 which resulted in the loss of one of the four potential N-linked 
glycosylation sites identified in case of the SIB fragments of the other BCV 
isolates. Residue substitutions lie for Ser for BCQ-20 and BCQ-9, and Thr for 
He for BCQ-9 also resulted in increased hydrophilicity (Fig. 3) and surface 
probability in this region. 

The phylogenetic tree in Fig. 4 shows the relationship of the four Quebec 
clinical BCV isolates to the published homologous sequences from BCV-F15 
[3], BCV-MEB [1], BCV-LY, BCV-L9, BCV-VACC [37] based on their SIB 
fragment genes. The homologous sequences of BCV-QUEB [20], a cell culture- 
adapted enteropathogenic strain isolated in Quebec in 1977, was also included 
in this comparison. Phylogenetic analyses by Unweighted Pair Group Method 
with Arithmetic Mean (GeneWorks) classify recent Quebec clinical isolates in 
a distinct sublineage. 

Comparison of the proteolytic cleavage site sequences 

The third region that was examined contains the coding sequence for the 
proteolytic cleavage site. The prototype BCV cleavage site sequence consists 
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of six residues KRRSRR beginning with aa 763, and predicts a cleavage between 
aa 768 and 769. Cleavage at the Arg-Ala bond by a host cell trypsin-like enzyme 
generates the SI and S2 subunits of the S protein [1]. The sequence of the 
predicted proteolytic cleavage site (KRRSRR) was conserved among all the 
isolates studied. However, three of the Quebec clinical BCV isolates examined 
had modifications in aa proximal to the cleavage site. The isolates BCQ-2070, 
BCQ-20 and BCQ-9 had a Ala-769 to Ser change immediately upstream the 
cleavage site. 


Discussion 

Considering the complexity of the pathogenicity and tropism of BCV and other 
coronaviruses, such as turkey coronavirus or human coronavirus HCV-OC43, 
comparative analysis of different clinical isolates will probably not produce 
meaningful results. However, with accumulating data on genomic and antigenic 
characterizations, immunodominant and functional regions can be selected for 
sequence analysis. In the present study, we chose four Quebec BCV isolates 
recovered from cases of epidemic diarrhea in neonatal calves in Quebec dairy 
herds between 1987-1990 to analyze their sequences in a selected genomic area 
which codes for immunologically and functionally important domains. Nucleotide 
sequence analysis of the four clinical isolates indicated that their SIB gene 
fragments were highly conserved. When compared with the prototype BCV- 
MEB strain, the SIB gene fragments of these Quebec isolates differed by only 
2.5 to 3.5%, the highest concentration of nt divergence being located in the 
polymorphic region. This region has greater than 94% nt identity with that of 
the prototype Mebus strain, and apparently, the differing nt were not distributed 
randomly over the sequences. The mutations have different effects on the 
encoded aa sequences. Over the SIB sequence most of mutations were silent, 
in contrast, in the polymorphic region most of the nt substitutions resulted in 
changes of aa residues. The SIB deduced aa sequences were more than 97% 
conserved among the strains, however, the polymorphic region was only 91% 
conserved. Recently, comparisons of the S gene sequences of six highly virulent 
(including BCV-F15) or avirulent (including BCV-MEB, and BCV-QUEB) 
strains have revealed more than 98% sequence identity in paired comparison, 
and most of the substitutions occurred in the SI subunit and were distributed 
randomly over the entire sequence [37]. The polymorphic regions were almost 
completely conserved among all strains, except for the virulent BCV-F15 strain 
[37]. In the present study, a closer genomic relatedness was demonstrated 
between BCV-F15 strain and the four Quebec clinical BCV isolates. At least 
twelve aa substitutions have been identified between the virulent and avirulent 
(BCV-MEB and BCV-QUEB) groups in the highly polymorphic region of the 
SI peptide (Fig. 2), suggesting that aa changes in this particular region of the 
S glycoprotein may be related to BCV virulence. However, the results obtained 
do not exclude that other-regions of the BCV genome may be also involved 
in the virulence or tissue tropism. 
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It has been previously demonstrated that antigenic domain II, identified 
within aa residues 517 to 720, is functionally and immunologically important 
[36], The epitopes that mapped in this domain are capable of inducing 
neutralizing antibodies [34, 36] and, therefore, represent targets of immuno¬ 
logical pressure. This antigenic domain, together with the polymorphic region, 
demonstrated a variability among the Quebec clinical isolates. Although the 
number of aa substitutions was relatively small, the fact that most substitutions 
appeared to cluster in an hypervariable region of the S gene, suggests these 
differences may be significant. The present data are in agreement with recent 
findings using various immunological assays that have shown that the four 
BCV isolates tested, reacted differently with neutralizing MAbs directed against 
the SI subunit of the Mebus strain [19]. 

Three of the Quebec BCV isolates examined could be differentiated from 
the prototype BCV-MEB strain by the nt sequences proximal to the proteolytic 
cleavage site. These isolates had a Ala to Ser change at residue 769 immediately 
upstream the cleavage site. This change maintained the number and pairing of 
basic residues of the prototype sequence. Changes in amino acids proximal to 
the cleavage site of the HA molecule of the influenza A virus have been 
demonstrated to modulate its cleavability [23]. In addition, changes from Phe 
to Ser immediately following the proteolytic cleavage site of the F protein of 
Sendai virus conferred the acquisition of susceptibility to plasmin cleavage; this 
suggests that regions of the F protein contiguous to the cleavage site can also 
affect its cleavability by interfering with the efficacy of host proteases involved 
in the maturation of the viral envelopped glycoprotein [14], However, in the 
present study, the Ala to Ser change immediately following the cleavage site 
of the S glycoprotein of three of the Quebec BCV isolates did not appear to 
modulate cleavability, and no correlation was found between this substitution 
and the rate of viral replication and type of cytopathic changes induced in 
HRT18 cells (data not shown). 

The genomic variations which have been identified in the present study 
seemed to represent characteristic features of the Quebec BCV isolates. Con¬ 
sidering the fact that these isolates were not passaged more than five times in 
HRT-18, and because of the finding of limited identity between Quebec BCV 
isolates, the possibility that our sequence data contain only aa substitutions 
arising after cell culture passages may be excluded. However, such aa substitu¬ 
tions may have occurred in cases of BCV-MEB and BCV-QUEB isolates for 
which numerous passages were done in cell lines other than HRT-18 and from 
heterologous animal species [7, 27]. The human rectal adenocarcinoma cell 
line has been reported to be susceptible and permissive to human, canine and 
turkey enteric coronaviruses [6, 8,16,17]. This human cell line grows in culture 
in tightly packed colonies of cells that have retain some characteristics of mature 
enterocytes [32]. The adaptation of field strains of enteric coronaviruses to this 
cell line usually requires only few passages. For RNA viruses, it is expected 
that numerous cell culture passages can result in genomic variations which 
may occur as point mutations randomly distributed throughout the entire virus 
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genome, rather than clustered in a specific genomic region. In the present study, 
the sequence heterogeneity was located mostly in the polymorphic region of 
the SI glycoprotein. 

In summary, our data have shown that Quebec clinical BCV isolates are 
genetically divergent from the reference strains [1, 3, 20, 37], and that the 
polymorphic region demonstrates a high variability. The antigenic as well as 
genetic relationships among the currently circulating strains of BCV provide 
useful information for determining the most appropriate composition for the 
coronaviruses vaccines. Continued surveillance of these viruses will be necessary 
to detect the emergence of new lineages of BCV with altered antigenic properties 
and assure that the future vaccines contain the most appropriate strains of BCV. 
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